Quantum Dynamics and Statistics of 
Atom-Molecule Bose-Einstein Condensate 
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Based on a two-mode boson model, we study nonclassical properties of the atom-molecule Bose- 
Einstein condensate. The effects of nonlinear collisions on the dynamics of the molecular formation 
is studied both in classical and quantum treatments. We find that the conversion from atoms 
to molecules can be suppressed strongly due to nonlinearity-induced localization of the atomic 
population. In addition, we study statistical properties of the atom-molecule condensed system by 
calculating the intensity correlation functions numerically. We find that the effect of nonlinearity 
leads to the appearance of superchaotic molecular pulses, while maintaining the atomic field sub- 
Poissonian. The joint quantum statistical properties of the atoms and the molecules always show 
anti-bunching. 
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Recently, possibility of preparing atom-molecule Bose- 
Einstein condensate (AMBEC) has attracted wide atten- 
tions [1-4]. Condensed bosonic atoms can be converted 
to a molecular condensate by using either the photoasso- 
ciation process [1,5,6] or the so-called Feshbach resonance 
method [2,3,7,8]. Wieman's group in JILA [7,8] measured 
the first-order temporal mutual coherence of atoms and 
molecules by using two-photon Ramsey experiment. Sim- 
ilarly with previous works [9,10], almost coherent Rabi 
oscillations (Ramsey interference pattern) were demon- 
strated. Theoretically, their experimental results were 
reproduced partially by using coupled Gross-Pitaevski 
(GP) equations [11-14]. According to the mean-field the- 
ory (MFT), large-amplitude coherent oscillation between 
the two-field modes is expected. Vardi et al. [15] found 
that quantum-field solutions modify the Bose-enhanced 
factor of the oscillation frequency. Moreover, many-body 
quantum effects lead to the appearance of collapse and 
revival of the coherent oscillations [5,15]. 

In this letter, we study nonclassical properties of the 
AMBEC based on a two-mode bosonic model [5,15]. The 
role of nonlinearity on quantum dynamics of the AMBEC 
is studied. We find that nonzero interspecies and in- 
traspecies interactions result in modulational instability, 
and modify the dynamics of the AMBEC. The conversion 
of atoms into molecules is shown to be dramatically sup- 
pressed due to nonlinearity-induced localization of the 
atomic population. We also study the statistical prop- 
erties of the AMBEC by calculating the intensity corre- 
lation functions numerically. Our results show that due 
to the nonlinearity, the initial sub-Poissonian molecular 
field is transformed into superchaotic molecular pulses, 
whereas the atomic field remains in the sub-Poissonian 
region. In addition, the joint quantum statistical proper- 
ties of the atoms and the molecules always exhibit anti- 
bunching. 

We consider the atom-molecule condensate system 
coupled via one-color photoassociation or Feshbach reso- 



nance. The total system can be described phenomenolog- 
ically by a two-mode bosonic Hamiltonian [5,15] (h = 1): 
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where 5 denotes the detuning between the molecular and 
atomic field. The atomic- and molecular-field operators, 
a and b obey the standard bosonic commutation rela- 
tionship. The atom-atom, molecule-molecule, and atom- 
molecule elastic interactions are described by the s-wave 
scattering strengthes A a , At, and X a b (Ab a ) [2,6], respec- 
tively. It should be mentioned that Af, and X a b are still 
unknown [4,6]. Following Refs. [16,17], we take A& ~ 2A Q 
(~ 10~ 3 g) and X a b ~ — 1.5A a . Due to the conserved 
particle number N — n a + 2fib, where n a — a*a and 
hb = b'b, the Hilbert space is spanned by the basis [18]: 
\(t>n) = \N - 2n) a \n) b for 71 = 0, 1,..., [N/2], where [z] 
denotes taking the greatest integer of any real number 
z. The dynamics of the AMBEC for arbitrary initial 
states can be solved by using the Runge-Kutta numeri- 
cal scheme. 

Before investigating the quantum dynamics of the AM- 
BEC, we first derive the classical equations of motion and 
discuss their steady states. The GP-type equations for 
the atomic and molecular fields can be obtained from the 
Heisenberg equations by replacing a — > (d) = \J {N)<fr a 
and b — > (b) = \J (N)$>b, where the wave functions for the 



two-field modes $i 



U 



b) obey the nor- 



malized condition: |$ | + 2|<5>{,| 2 = 1. To get the steady 
states and the modulational instability conditions [19], 
we introduce two conjugate variables [20,21]: x = |$ | 
and if — 2<j) a — 4>b. They obey the canonical relations: 
dx/dr = —dH/dtp, dp/dr = dH/dx, and the canonical 
Hamiltonian is given by 

H = Ax — ^-x 2 — x\f\ — xcos(ip), (2) 
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where the scaled time r = \f2Gt. In derivation of the 
above equations, we have introduced 

A = (V2G)- 1 ((5-A ab + A h /2), (3) 
A=(G/V2)- 1 (A a -A ab + A () /4), (4) 

where G = \J (N)g and A^ = (N)Xij. The stationary 
solutions obey (dx/dr) Xo = [dip/dr)^ = 0. From the 
relation x n ^l — xo sin(</?o) = 0, we find that there are 
three kinds of fixed points: (i) pure molecular phase with 
xq = 0, (ii) in-phase type of steady state with ipo = 0, 
(iii) out-of-phase type of steady state with ip = n [16]. 
The value of x determined by a cubic equation 

A - Ax - 1 _^/, 2 cos(yo) = (5) 

can be solved numerically. It was shown that the atomic 
phase at |A| — > oo can be converted to the molecular 
phase or vice versa by varying the detuning A adia- 
batically [5,21]. In the intermediate region, there are 
two-type atom-molecule coexisting phases: the in-phase 
steady state for ip = in the parameter region A < 1, 
and the out-of-phase one for ip = ir in A > — 1 region. 

Stability analysis of the steady states can be performed 
by introducing small perturbations: x = xq + 5x and 
f = fa + 8<P- A coupled linearized equation for the per- 
turbations gives the eigenfrequencies 

oj 2 = ojq — Ax \/1 — xo cos(ipo), (6) 

where lu = \J x (l — 3x /4) cos 2 (<po)/(l — x ). Fot a 
negligible nonlinear interaction A = 0, aroused from ei- 
ther Aj = Aij = or A ab = A a + A b /4, both of the 
two steady states (920 = 0, ir) are stable [21]. However, 
for nonzero A, the above two steady states can be mod- 
ulationally unstable when the eigenfrequencies become 
imaginary [22], i.e., uo 2 < 0. In this case, the perturba- 
tions in the AMBEC undergo exponential growth with 
the growth rate being the imaginary part of u>. We find 
that, for (po = 0, the MI occurs in the positive A region 
with A > 1 and A < 1. For tp = n, the AMBEC is un- 
stable in the negative A region with A < — 1 and A > — 1. 
The MI region (and also the dynamics) of the AMBEC is 
invariant under the transformation A — > —A, ip — ► it — ip, 
and A — > —A. Thus, in the following we only consider 
the case of positive A. 

Dynamical evolution of the AMBEC can be described 
by the fraction of atoms converted to molecules [5,19], 
/m = 2(6^6)/(A r ), which measures the ratio of the atom 
number in the bound state to the total atom number. 
In the MFT, f M = 2|$ fc | 2 . For the negligible nonlinear 
interaction A = 0, there are analytical solutions for Jm 
[23]. Considering the initial condition x(Q) = 1, then 
/m = 2/isn 2 (v / ^r/v / 2, m), where sn(^, m) is a Jacobi 
elliptic function with the elliptic modulus m = /1//2, 
and /i, 2 = I [A 2 + 2 =F ^A 2 (A 2 +4)]. In particular, for 



the exact "resonant" case with A = 0, /m = tanh 2 (r/2) 
since m = 1. Starting from the zero initial value, Jm 
grows quickly to the maximum value 1 without any oscil- 
lation (the dotted line in Fig. 1 (a)), which corresponds 
to a complete and irreversible molecule formation pro- 
cess. Nonzero A, however, results in periodic oscillations 
of Jm with smaller amplitude. If A > A c = 1/V2, the 
amplitude of Jm is always less than 1/2. 

To investigate the effect of nonzero A, we need to 
know the values of the two-body interaction strcngthes. 
However A b and X a b have not been determined even in 
low-energy region [4,6]. Following Ref. [17], we adopt 
A b = 2\ a and A ab = — 1.5A a . Nonzero A can be achieved 
by tuning A , such as A a = 1.0 x 10~ 3 .g (A = 0.4243), 
A Q = 2.22 x l0- 3 g (A = 0.9419), and A a = 4.0 x lO^g 
(A = 1.6971). All the above parameters are in the sta- 
ble region. The MFT results show that the effect of 
nonlinearity leads to periodic oscillations between atoms 
and molecules. As shown in Fig. 1 (a), the ampli- 
tude of /m decreases monotonically with the increase of 
A. If A > 0.94, the whole curve of f M is below 1/2, 
i.e., no more than 50 percent of atoms are converted 
to molecules. This nonlinearity-induced localization of 
the atomic population is similar to the macroscopic self- 
trapping in the two-component BEC [20]. However, un- 
like the two-mode BEC system with Josephson-like cou- 
pling, the oscillation amplitude in AMBEC reduces con- 
tinuously from unity to zero. 

In Fig. 1 (b), we evaluate numerically the quantum 
solutions of Jm for an initial number state |-0(O)) = 
\N) a \Q)b- In the absence of nonlinearity A — and 
A = 0, the quantum solution of /m (the dotted line 
in Fig. 1 (b)) breaks away from the MFT result around 
the point r ~ 7.03, where all the particles arc in molec- 
ular state. This departure can not be explained in the 
framework of the MFT. A semiclassical analysis of the 
AMBEC in the vicinity of the all-molecules phase shows 
that there are two contributions to the atomic occupation 
[15]: (i) spontaneous dissociation into atomic vacuum, 
which is proportional to \/~N sinh 2 (r/2), (ii) stimulated 
part proportional to (n a ), which is negligible near the all- 
molecules phase. The Bose-enhanced spontaneous disso- 
ciation into the atomic vacuum results in collapse and 
revival of the atom- molecular oscillations [5,15], which is 
just similar with the phenomenon observed in the two- 
photon down-conversion process [18]. 

We further consider quantum solution of /m for the 
case of A 7^ 0. As shown in Fig. 1 (b), we find that the 
quantum results almost coincide with the corresponding 
classical solutions shown in Fig. 1 (a). Both of them 
exhibit "periodic" oscillations with decreased amplitude. 
The agreement of two solutions originates from the fact 
that the trajectories do not reach the all-molecules phase, 
and the contribution of quantum fluctuation of atomic 
mode is not significant compared to the coherent stimu- 
lated process. In Fig. 1 (c), we compare the quantum 
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solution (red solid-line) with the MFT solution (dotted 
line) for A = 0.4243. We find that the quantum solu- 
tion is in fact the collapsed oscillations. From Fig. 1(b), 
we also find that the collapse time is enhanced signifi- 
cantly by the nonlincarity, i.e., the coherent oscillations 
are damped more slowly with the increase of A. 



then the field is in the so-called sub-Poissonian. The 
fields in a sub-Poisson distribution always exhibit non- 

(2) 

classical anti-bunching phenomena. <?j — 1 represents a 

(2) 

coherent state with Poisson distribution, while g) > 1 
is characteristic for a super-Poisson distribution. In par- 

(2) 

ticular, for thermal or chaotic field, = 2. 
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FIG. 1. (Color online), /m versus the scaled time r for (a) 
mean-field solutions, (b) quantum solutions, and (c) compar- 
ison of the two predictions for A = 0.4243. The parameters 
are A = 0, N = 10 4 , and the nonlinearity in (a) and (b) 
are: A = 0.0 (the dotted lines), A = 0.4243 (the red lines), 
A = 0.9419 (the blue lines), and A = 1.6971 (the black solid 
lines). 
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FIG. 2. (Color online), and g^ versus the scaled 

time t for various A. From top to bottom: (a) A = 0, (b) 
A = 0.4243, (c) A = 0.9419, and (d) A = 1.6971. The red 

(2) (2) 

solid lines are for ga and the blue dotted lines for <j£ . Other 
parameters are the same with Fig. 2. 



The observation of the first-order temporal coherence 
between atoms and molecules opens up the possibility 
to form molecular condensate. However, there remains 
many open questions to be answered [7] , such as the de- 
tailed quantum state of the molecules, and their higher- 
order coherence functions in the temporal domain. Re- 
cently, Meiser et al. [24] studied the short-time dynamics 
of molecule formation and the second-order correlation 
function by using perturbation theory. They found that 
the second-order degree of coherence of the molecular 
field for an initial number state approaches unity in the 
limit of large N. Motivated by their work, we calculate 
the equal-time intensity correlation functions beyond the 
short-time limit: 



It is known from quantum optics that if < g < 1 



Starting from the number state, the initial stage of 
the second-order coherence functions of two field modes 
obey g ( a ] = I - l/N and gf ] = 1 - 4/N + (D(N~ 2 ) 
[24]. Both of them are less than 1. For zero nonlin- 
earity A = 0, as shown in Fig. 2(a), we find that the 
two field modes are transformed from sub-Poissonian to 
super-Poissonian, successively. Similarly with previous 
discussions, the AMBEC system with zero nonlinearity 
shares the same model of degenerate parametric amplifier 
[25]. Quantum fluctuations of the atomic vacuum leads 
to the appearance of super-Poissonian, and the peaks 

(2) 

of <7a near the points r ~ 7.03 and 21.0, where the 
particles are almost in the molecular mode. The mag- 

(2) 

nitude of g a > 2 indicates that the atoms generated 
from the photodissociation exhibits superchaotic (strong 
bunching) behavior [26]. 

(2) 

The second-order coherence function g - for nonzero A 
are calculated in Fig. 2 (b) , (c) and (d) . We find that the 
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(2) 

effects of nonlinearity result in sharp peaks in g b with 
their positions at the minima of /m- Their magnitudes 
are larger than 2 but not divergent due to the nonzero 
minimums of the mean molecular number. Our results 

(2) 

imply that the initial sub-Poissonian statistics of g b is 
transformed into a super-Poissonian, and thereby a su- 
perchaotic molecular pulse can be realized at certain val- 
ues of time. The second-order coherence functions of the 
atomic field are also plotted in Fig. 2 (the red lines). We 

(2) 

find that super-Poissonian statistics of g a is fully sup- 
pressed for A ^ 0. More specifically, the intensity corre- 
lation function of the atomic field oscillates between sub- 
and super-Poissonian with very small amplitude. Fur- 
ther increase of A, such as to 1.6971, < 1, which 
implies that the atomic field is always in sub-Poissonian. 
The reduction of the super-Poissonian of atomic field to 
sub-Poissonian comes also from the suppression of spon- 
taneous atomic emission. The increase of A dominates 
the contribution of the coherent stimulated process. 

Finally, we study the joint quantum statistical proper- 
ties of the atoms and the molecules, defined by 



o (2) 

9ab 



We find that g ab (r) is always less than 1. Physically, 
such a result means that the atoms and the molecules 
do not tend to be created simultaneously, i.e., the anti- 
bunching between the two field modes. Similar results 
have been addressed for optical field in the two-photon 
down-conversion process [18]. 

In summary, we have studied the role of nonlinearity 
on quantum dynamics and statistical properties of the 
AMBEC. The former problem is related to the atom- 
molecule conversion rate. We have shown that even in 
the simplest two-mode model without "rogue photodisso- 
ciation" [27], the conversion can be strongly settled down 
due to the nonlinear analog of macroscopic self-trapping 
effect. We also studies statistical properties of the AM- 
BEC by calculating the intensity correlation functions 
numerically. Our results show that the nonlinearity A 

(2) (2) 

in the stable region can modify both g a and g b sig- 
nificantly. The initial sub-Poissonian molecular field is 
transformed into superchaotic molecular pulses, whereas 
the atomic field remains in sub-Poissonian region. Our 
result also show that the joint quantum statistical prop- 
erties of the atoms and the molecules always exhibit anti- 
bunching. 

It should be mentioned that Prataviera et al. [28] have 
studied higher-order mutual coherence of two chemically 
different species: the light field and the matter field. 
Their scheme make use of the combination of the tradi- 
tional absorption photodetector and a matter-wave de- 
tector relied on photoionization of the atoms. Simi- 
larly, photoionization of the diatomic molecules can be 



(2) 
9ba 



{n a ) {fib)' 



(8) 



proposed in the atom-molecule photoassociation exper- 
iments to measure the molecules and their statistical 
properties directly. 
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